String solutions in Chern-Simons-Higgs model coupled to an axion 
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We study a d = 2 + 1 dimensional Chern-Simons gauge theory coupled to a Higgs scalar and 
an axion field, finding the form of the potential that allows the existence of selfdual equations and 
the corresponding Bogomolny bound for the energy of static configurations. We show that the 
same conditions allow for the N = 2 supersymmetric extension of the model, reobtaining the BPS 
equations from the supersymmetry requirement. Explicit electrically charged vortex-like solutions 
to these equations are presented. 



I. INTRODUCTION 



There has been recently a revival in the interest on Abrikosov-Nielsen-Olesen flux tubes in connection with su- 
' pergravity models and the realization of cosmic superstrings. Indeed, although such cosmic size superstrings were 
originally excluded in the context of perturbative string theory 0| , it became clear recently that the question should be 
revisited as the relevance of branes and new kind of extended objects was understood (see 2] and references therein). 
In order to understand the nature and structure of such stringy cosmic strings, the embedding of BPS objects in 
supersymmetry and supergravity models has become an active area of research so that the properties of BPS solitons, 
their connections with the supersymmetry algebra and their cosmological applications have been discussed by many 
■ authors [^-[T^. 

Having in mind the study of BPS solitons in a string theory context, it is natural to consider models where an 
axion field is included. In particular, an = 1 globally supersymmetric model in d = 3 + 1 dimensions consisting 
of an axion superfield S coupled to WaW", with Wa the chiral superfield strength, was analyzed in 12] and finite 
energy cosmic string solutions were constructed. Also, the impact of axions on dynamics of a d = 3 + 1 Yang-Mills 
theory supporting non-Abelian strings has been analyzed in [iq . 

It is the purpose of this work to consider similar issues in d = 2 -I- 1 space-time dimensions for which a rich variety of 
flux tube solutions exists already when the axion field is absent. Indeed, when gauge fields with dynamics governed by 
a Chern-Simons action are coupled to charged scalars with an appropriate sixth order symmetry breakin g po tential, 
the model admits BPS equations with vortex- like solutions carrying both magnetic flux and electric charge |17|-[T^. It 
should be stressed that in the absence of the Chern-Simons term, electrically charged vortices with finite energy (per 
unit length) do not exist jl9|- Hence, the model we are interested in, could show novel aspects of charged string like 
configurations when an axion is present, in particular with respect to their application to cosmological problems. This 
in the perspective that at high temperatures, a relativistic four dimensional quantum field theory becomes effectively 
three dimensional. 

The paper is organized as follows: we start by considering, in section 2 a purely bosonic model with couplings and 
potentials chosen so as to have a gauge invariant action leading to non-trivial Bogomolny equations. These equations 
are the natural extension, when an axion field is present, of those found in ^3|"0I- We present in section 3 the 
supersymmetric extension of the model and we establish the connection between supersymmetry and BPS equations 
when the axion field is included. Numerical solution of such first order coupled equations are presented in section 
4, where we discuss, in particular, novelties in the vortex solutions resulting from the presence of the axion field. A 
summary and discussion of our results are presented in section 5. 
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II. AXION COUPLED TO A CHERN-SIMONS-HIGGS SYSTEM 

The coupling of an axion to a gauge field with dynamics governed by a Chern-Simons action poses some problems 
[l^-ll^l- To discuss how can they be overcome, let us consider the following (2 + l)-dimensional bosonic action, 



S 



d'x {^/(s)ImD^5F^ + \D^^\^ + K"is)\D^S\^ - W{<P, S)] (1) 



where (j) and S — s + ia are complex fields, is a U{1) gauge field, and F^^ is defined as, 

= e'""'d^A^ (2) 

W{(t), (j)* , S, S*) is a potential term, / and K are arbitrary functions of s, the real part of S, primes stand for d/ds, 
K is a constant and 5 is a dimensionless parameter (which in the 4 dimensional case is related to the Fayet-Iliopoulos 
term and the Planck mass). Finally, is the covariant derivative acting on the fields 4> and S according to 

D^S - d^S + 2i5A^ (3) 

As done in for the 3 + 1 model, we shall identify a in J^l with the axion field and s with a dilaton field. Note that 
because of the definition of the axion covariant derivative, the Chern-Simons term appears in action multiplied by 
the factor /(s), 

Scs[A,s] = ^J d'xf{s)e^'"'A^d,A, (4) 

The action |^ is invariant under the local transformation 

5 ^ S-2iSA{x) 
S* S* + 2i6A{x) 

^ e'<=^(^)cj) (5) 
A^ ^ Af, +df,A{x) 

The time component of the gauge field equation of motion is the Chern-Simons version of the Gauss law and can 
be used to solve for giving 

K B , , 

" " 2(e2|0|2 + 4(52if") ^ ' 



where 



B = fis)F,y - e^^{A, + la,a)9,/(s) (7) 



The energy can be found from the energy-momentum tensor obtained by varying the action with respect to the metric, 

5S=\j Sx^T>^^5g^, (8) 
Integration of the time-time component T'^" gives 
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After some work, this expression can be written in the form 

K"\D.^S±iDyS\^ 



£ ^ ld^x\ \D.:,(t>±iDy(t>\'^ 



4(e2|,/)|2 +4(52i^") 



± e( 



9jlog/(s) 



(10) 



We thus see that the first three terms in ea. pU|) have been accommodated as perfect squares. This, together with an 
appropriate choice of the potential W so as to cancel the fourth and fifth terms, would lead to a Bogomolny bound 
for the energy given by the magnetic flux $ appearing in the last term, 



or 



$ = l-nn 



(11) 



(12) 



for gauge fields with topological number n e Z. There is however the sixth term in (|10|l with no definite sign preventing 
the obtention of a bound. Only if we put /(s) = 1, which corresponds to a normal Chern-Simons action for the gauge 
field (see eq.(@J)), this term vanishes. In that case one does have a bound. 



f > ±e|0o|'$ = 2^e|(/)oP|n| 



whenever the potential is chosen as 



(13) 



(14) 



The bound is saturated by fields obeying the self-duality equations 

^'a;^ = TiDycj) 
D,S = TiDyS 

K^F.y = t2 (e'|</.p + AS'K") (e(|0|2 - j^oH - ^SK') (15) 

where the upper (lower) sign corresponds to positive (negative) values of (f>. 

As it is well-known, the presence of a Chern-Simons term forces a relation between magnetic flux and electric charge 
|22|. This makes the Chern-Simons vortices both magnetically and electrically charged. To see this phenomenon in 
the present case let us write the gauge field equation of motion, for the case f{s) = 1 in the form 



(16) 



with = (p, J) the conserved matter current and p the electric charge density, 

= -2 {e^\(l)\^ + 4K"6^) - ie{(f>d''(f>* + (/)*a^0) - ASd^^a (17) 
We then see that eq.® can be rewritten in the form 

P = -kB (18) 
so that the usual relation between electric charge Q = J p d^x and magnetic flux in Chern-Simons theories holds, 

Q = (19) 
Note that both the Higgs scalar and the axion contribute to the electric charge. 
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III. SUPERSYMMETRIC EXTENSION 



The SUSY extension of the Chcrn-Simons-Higgs system with a sixth order symmetry breaking potential was ana- 
lyzed in |23|. Let us study now the the case in which the axion field is also present. 
Consider the following d = 3 action, written in terms of superfields as 

5susY ^ -IJ d'2;d20(^F(I] + Et)im(^V„.I]) M^" 

+ jy($t, $)(v'^$)tv,$ + K,si^, s^)|VaS|' 

+ 2T/($t,$,S,Et)) (20) 
The index o takes the values 1, 2 and r^, $ and E are spinor, complex scalar and axionic superfields. Also 

Va$ = {Da ~ ieTa)<i> 

where /i = 0, 1, 2. 

The lowest component of F", $ and E are respectively, the gauge field, the Higgs field and S = s + ia with a 
the axion field. Concerning F, H, K and V, they are functionals of superfields to be fixed later. Subindices in these 
functionals indicate derivatives, thus Ksg — dsds*K = c^s^gi^le^o and so on. These functionals should be chosen so 
that the supersymmetric action H2(J|I is invariant under the supergauge transformations, 

Ta ^ r, - DaA 

E -> T,a-25A (21) 
for any real scalar superfield A. Written in the Wess-Zumino gauge, i.e., 

ra|e=o = D'^Ta\e=o = 0, 

the spinor superfield Ta is given by 

where A(a;) is a Majorana spinor, the photino. Then, the spinor field strength, defined as 



1 



Wa - -D''DaTt 



takes, in terms of component fields, the form 



and satisfies the Bianchi identity. 



D'^Wa = 0. 



The complex scalar superfield is defined as 

$(a;, 6) = (f>{x) + e^ijaix) - e^F{x) 

where cj) stands for the Higgs complex scalar field, ipa is a Dirac bispinor, the higgsino and -F is a complex auxiliary 
field. Finally, the superfield E which contains the axion S as its lowest component, can be reduced to a complex 
scalar field O by exponentiation. 



e = 
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The supersymmetric transformations H21() take out from the Wess-Zumino gauge. One can however implement a 
composition of SUSY and gauge transformations such that the Wess-Zumino gauge remains vahd. To do that, the 
new SUSY transformation for scalar and spinorial superfields are, respectively, 

S^^Ta = Iv'QbTa + Dak 

,5^2$ = irfQt,'^ + ik<^ (22) 
where Qa = ida + 6''' (7^) 6a 9^, is a Majorana spinor and the real scalar superfield K is defined as, 

k = ie''{-f'^)abv''A^+0^X^7Ja. 

Let us restrict the action (|20|) to the case in which 

+ =i7($t^$) = 1 (23) 
Then, written in components, the action takes the simple form 

SsusY — Sb + Sf (24) 

with the pure bosonic action given by 

Sb = J d?x l^e'^-^^F,, + |D^0|2 + R\b^S\^ - \V^\'' - ^\Vs?^ (25) 
while the fermionic one is given by 



+ ^p^p ( -101^ -v^^\+xx[ —R' ~ V,, 



1 , (R, A 1 452 , 

+ Vi*x [^rR ~ v^s^ + v^x* - ^0.) } , (26) 

where we have redefined R = Kss- Here we have taken into account that 

I = j d^xffe [D^Y.) Wa 
is a surface term which does not modify the equations of motion 

/ = j d^xd^OD" {T,Wa) = J d^xD'^D'' (SW^a)le=o 
= z j d'xd^D' {YWa)\o^o . 

and can hence be neglected. 

In order to extend the A'^ = 1 supersymmetry to an = 2 supersymmetry, we shall follow 23] and allow the 
transformation parameter 77 to be complex, i.e., rf ^ ry. Now, this is the same as making an = 1 SUSY transfor- 
mation followed by a U{\) fermion phase rotation. Thus, the new transformation for fermions will be tpa e'^V'a 
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and ip* e~*"V'a ^-^^ the same is valid for Xa and Xa- The new SUSY transformations act then as rotations on the 
fermions and one can then see that the only terms which do not respect the extended SUSY invariance are those on 
the last three lines in H2f)|l. Hence, in order to get an = 2 supersymmetric model we need that, 



4^2 



-R' - Vss - 



p2 


= 0, 


2Se 


K 




K 


4^2 




26e 


i?2 


= 0, 




K 


K 



K 

where V — V{u^ v), with u — (f)*(j) and v — S + S* . These equations imply that, 

e 



Vu 



K 



{eh{u) - 25r{v))H{u) 



where 



dv 



25 



{26r{v) - eh{u))R{v) 



r{v) = R{v) 



du 



h[u) 



(27) 



We obviously have h — u~ uq = |0|2 — |0o|2, and since S + S* ~ 2s, then R = K" and r = 2K' where primes stand 
for derivatives with respect to s. From H25|) the potential is, 



W 



4(e2|0|2+4<52i^") [e( 



\M^)-A5K'] 



which is exactly what we obtained in l|14(l . 

In order to get the Bogomolny bound and the self-dual equations one can analyze the supercharge algebra as in 
Alternatively, one can directly consider the component field SUSY transformations ((5^X = ry"(5aX), 



SaF^i{^'')lD,,ijl, + 2Xa(j), 



SaS = -Xa 
SaXb = <^abJ - i{Y)abD^S 

SaJ = i{lnl^^.Xb + 2Xa 



(28) 



and their complex conjugated {SX'' = ri*'^6aX'') and reobtain Bogomolny equations (|15|l just by putting all fermion 
fields to zero in 128|1 and then ask the SUSY transformations for ipa and Xa to vanish once the auxiliary fields have been 
written in terms of dynamical fields using their equations of motion. The first condition corresponds to a restriction 
to the bosonic sector, the second one implies that physical states are supersymmetry invariant. 



IV. VORTEX-LIKE SOLUTIONS 



We present in this section some vortex solutions to the BPS equations of motion. Wc choose for the Khaler potential 
the form 



= -M2log (5 + 5**) 



(29) 



As in we shall analyze separately two cases: first, we consider the case in which the vortex is supported by the 
Higgs field ("0-strings" solutions) , in the sense that at infinity it behaves as in the ordinary Nielsen-Olesen vortex, 
with its winding number linked to the magnetic flux. Then, we shall consider the case in which the vortex is supported 
by the s field, a solution that we shall call an "s-string" . In this case it is the axion winding number which is related 
to the magnetic flux. 

In the first case, in order to obtain 0-string solutions we make the ansatz |l2| 



(l>i = fir)e 



S = s{r) — 2iSm9 



v{r) 



(30) 
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where n is the topological charge of the Higgs and m is the topological charge of the axion. 
It is convenient to work with dimensionless variables by defining 



ar, 



a 



With this convention the equations read, 



x{t) = v{t /a) 
uir) = eS^^ s{t / a) 



\n\ yO ' ^ 

2 

- — (|m| — |7i| a;) 



1 + — 



4/3 

y 



T 

where f3 = M^/0q. From the first two equations we can integrate y(r) in terms of z{t), obtaining 

2/(t) = 2(|n| - |m|) logr - 21ogz(T) + k 



(31) 

(32) 

(33) 
(34) 

(35) 



with k an arbitrary integration constant. Thus, we end with a system of two first-order coupled differential equations 
for x{t) and z{t). 

The boundary conditions for equations H32() and (|34|l can be determined as follows. By equation H30|) . the function 
x(r) must vanish at the origin, so eq. H34|l implies that z(r) also vanishes at the origin as r'"'. Indeed, eqs. (|32|l and 
(Enj) imply that x ~ for m — and a; ~ r^/ log(T)^. 

For large r, the function x tends to 1, thus z also tends to 1 unless \n\ = \m\. In this last case, z zq, where zq 
is the solution of the algebraic equation 



4 + 



4/3 



k — 2 log Zq 



1 



(36) 



To solve the differential equations we employ a relaxation method for boundary values problem. Such a method 
determines the solution by starting with an initial guess and improving it iteratively. We present in figure Q some 
solution profiles for the gauge field x{t) and for the dilaton ?/(t). 



1 • 



0.8 



X{T) 




Figure 1: Plot of the gauge field x(r) and the dilaton j/(t) for n — 1 and m — (solid line) and for n = 2 and m — 1 (dashed 
line). In both cases /? = 1 and k — 10. 



Already from ansatz (|30|l one can see that the magnetic flux and electric charge ar quantized according to 

27r ^ 2ttk 

$ = — n , Q = n 

e e 



(37) 
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There is an interesting property, typical of Chern-Simons vortices that also holds in our model: both the magnetic 
field and the (radial) electric field are concentrated in rings surrounding the zeroes of the Higgs field. We show this 
behavior in figure El 
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Figure 2: Vortex magnetic and electric fields. The solid line corresponds to n = 1 and m = and the dashed line to 
n = 2, m — 1. In both cases /3 = 1 and k = 10. 

We also solved the BPS equations when the field s tends asymptotically to a constants ("s-string", In this 

case the magnetic charge is equal to the topological charge of the axion field, so we have as an appropriate ansatz, 

(j) = f{r)e"'^ S = s{r) - 2idme Ae = (38) 

With this ansatz the BPS equations take the form 

|to| V y J \ y J 

y' = -- |to| (1 -x) (40) 

T 

z' ^ -{\n\~-\m\x) (41) 

T 

Again, we can integrate 2/(r) in terms of z{t), obtaining 

2/(t) = 2{\n\ - |m|)logT - 21ogz(r) + k (42) 
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We see from this equation that for consistency, z(r) ~ ^n-m j^j. ^ qo in contrast with what happens for the 
^-string. Concerning the gauge field boundary condition, one has liuiT^ao x{t) = 1. We present in figure |31 some 
profile of the resulting s-string solutions. 



1.2 

1 

0.8 
0.6 
0.4 
0.2 



2 4 6 8 



Figure 3: Profile of an s-string solution. The solid line corresponds to the gauge field x{t), the dot-dash fine to the dilaton 
j/(r) and the dashed line to the Higgs field z{t). For this solution n = 1, m = 2, /S — 0.1, and k — 0.1. 



V. DISCUSSION 

We have been able to write first order Bogomolny equations for the Chern-Simons-Higgs-axion system provided 
the gauge field dynamics is governed by a normal Chern-Simons action and the potential takes the form (|f4|) . with 
coupling constants satisfying the usual relation. As in the 3 -I- f case discussed in the first order equation for 
the scalar S containing the axion field has the same form as the ordinary Bogomolny equation for the Higgs field. 
The sixth order symmetry breaking potential (necessary to accommodate the energy as a sum of perfect squares plus 
a topological term) is modified due to the presence of the axion. As it always happens when matter is coupled to 
a Chern-Simons term, magnetic flux and electric charge turn to be proportional (see ea. H19|) ). Were we considering 
non-Abelian gauge fields, then the proportionality constant k should be quantized, implying in turn quantization of 
both electric charge and magnetic flux at the classical level. 

The same constraints and results were found by constructing the SUSY extension of the purely bosonic model. In 
this case we have seen that the conditions on the potential and coupling constants arise when extending supersymmetry 
from iV = f to = 2 while the BPS equations can be inferred by asking physical states to be annihilated by SUSY 
generators. 

The numerical solutions to the BPS equations presented in section 4 show that the axionless string solutions found 
in are not much modified by the axion which, however, contributes to the electric charge of the string 

configuration in a clear way. These explicit solutions could be of relevance in the context of cosmic strings and, 
due to the coupling to the axion and their electric charge, their dynamics could be very different from that usually 
considered. This, and the extension to the non-Abelian case (where both the magnetic flux and the electric charge 
are quantized) are issues that we hope to address in a future work. 

VI. APPENDIX: CONVENTIONS 

We use the following representation of the 7 matrices, 

Thus, 




7 7 = 5 - «e 7a 
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The "metric" in spinor space is given by 

/ n 1 \ 

= (43) 



So, 



■/ 2n 1 \ 



In particular, 



We define the integration measure in the superspace, 

^ d^9 = -1 
Also, useful formulae are, 

Ba0h = ~Cah0^, B^d^ ^-C^G^, d°'B},=&lB^ = -B^B°- 
The 7 matrices are real and symmetric when lowering their indices, i.e.. 

Then, any vectorial representation can be written in terms of the 7 matrices with spinorial indices, i.e., 

Vab = iinabVf, 

where, "fab are imaginary and symmetric 2x2 matrices. Also, we can define the spacetime derivative in terms of 
spinorial indices, 

dab = il'')abd^ 

The spinorial derivative is define as, 

Then, the supersymmetric covariant derivative can be written in terms of this notation as, 

Da = da+ le^'dba 

Note that the derivatives and do not raise and lower indices in the same way as spinors, 

a. = (wA^c..(a„.^)A.a.A.c..5^ 



_dx^_d__ dx^_d_ _^_r' rid 
causing the appearance of a minus sign in the definitions of 3°^ and 9^ in (|45|l and (|46|l . 
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